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I.  INTRODUCTION 


Nonlinear  transient  analysis  of  structures,  especial ly  those  involving 
material  nonlinearities,  are  highly  complicated  even  for  the  uniaxial  case 
with  a  simple  constitutive  model.  In  the  case  of  the  beam  bending 
problem  for  instance,  there  is  a  single  component  of  strain 
along  the  longitudinal  beam  axis.  This  axial  strain  component  is  assumed 
to  vary  linearly  over  the  cross-section  by  virtue  of  the  assumption  that  plane 
sections  remain  plane.  However,  since  the  stresses  are  no  longer  proportional 
to  strains  in  the  inelastic  range  the  distribution  of  stresses  can  be  quite 
complex  over  the  cross-section  and  along  the  length  of  the  beam.  Thus,  for 
a  beam  which  has  yielded  in  a  certain  region  the  line  of  zero  stress  or  strain 
does  not  coincide  with  the  centroidal  axis  in  that  region.  Hence,  in  general 
the  centroidal  axis  is  not  an  axis  of  constant  strain  for  deformations  in  the 
inelastic  range. 

A  two-noded  beam  column  element  is  often  used  to  analyze  inelastic 
response  of  models  built  up  from  such  elements  in  conjunction  with  other 
elements.  A  linear  axial  displacement  field  and  a  cubic  transverse  dis¬ 
placement  field  are  used  in  arriving  at  the  stiffness  properties  of  such 
an  element.  If  the  centroidal  axis  is  used  as  the  reference  axis,  it 
is  not  possible  to  satisfy  equilibrium  in  the  inelastic  range.  The 
same  is  true  of  the  linear  elastic  range,  if  an  axis  other  than  the 
centroidal  axis  is  used  as  the  reference  axis.  Hence,  a  simple  linear 
elastic  analysis  with  an  arbitrary  reference  axis  can  be  used  to  demon¬ 
strate  the  point  under  consideration. 

Another  feature  of  the  nonlinear  analysis  by  the  finite  element  method 
is  the  necessity  of  integrating  complex  distributions  of  stresses  and 
strain  energy  densities  over  the  volume  of  the  element.  Because  of  the 


complexity  of  the  integrand  and  the  domain  of  integrations  recourse  to 
numerical  schemes  like  Gaussian  quadratures,  Newton-Cotes,  etc.  has  to  be 
made  in  order  to  obtain  approximate  estimates  of  stress  resultants  and 
total  energies.  The  number  of  integration  points  used  over  the  cross- 
section  and  over  the  length  of  the  beam  element  determines  the  degree  of 
approximation.  Accordingly  approximate  estimates  of  energies  yield 
solutions  of  varying  degrees  of  accuracy  which  depend  upon  the  order  of  the 
integration  scheme  used.  For  an  assess.ment  of  the  quality  of  such 
approximate  solutions  it  is  necessary  that  an  exact  solution  to  a  problem 
be  known.  Such  a  problem  along  with  its  accompanying  exact  solution  has  been 
outlined  and  a  rigorous  evaluation  of  the  sensitivity  of  the  quasi-static 
response  to  the  order  of  the  integration  scheme  is  made.  For  problems, 
especially  those  involving  transient  response,  for  which  no  exact  solutions 
can  be  obtained,  only  qualitative  estimates  of  sensitivity  can  be  obtained 
by  a  comparision  of  numerical  solutions  using  different  orders  of  integration 
schemes . 


II.  BEAM-COLUMN  ELEMENT 


In  this  section  the  two-noded  beam-column  element  of  Fig.  l.a  will  be 
shown  to  yield  inaccurate  results  when  used  to  model  a  beam  bending  elasti¬ 
cally  about  its  centroidal  axis,  but  analyzed  with  respect  to  another 
reference  axis.  With  the  finite  element  formulation  in  mind  the  present 
discussion  will  be  confined  to  Euler-Bernoulli  beams  subjected  to  concen¬ 
trated  shears  and  moments.  For  such  a  beam  the  total  potential  energy  of 
deformations  is  given  by 

^  =  If  I  2  dAdx  -  P,w(0)-P2w()i)-M^ 

0  A 

(1) 

where  u  and  w  are  the  axial  and  transverse  displacements  of  the  reference 
axis  and  z  is  measured  normal  to  the  .reference  axis.  Minimization  of  tr 
with  respect  to  u  and  w  yields 


El 


.  eza  .  0 

dx^  dx"* 
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I  =I-Z^A  (5) 

c 

I^  in  Eq.  (5)  is  the  moment  of  inertia  of  the  cross-section  about  the 
centroidal  axis  which  is  separated  from  the  reference  axis  by  Z. 

Equation  (3)  integrates  to 

u  =  z^+cix+c  (6) 
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Assunia  such  a  baam  1s  to  be  modeled  using  the  beam-column  element  of 
Fig.  l.a.  This  finite  element  assumes  a  cubic  transverse  displaceinent 
of  the  reference  axis.  For  satisfaction  of  equilibrium.  Eq.  (6) 
requires  that  the  axial  displacemant  vary  quadratical ly  over  the  length 
of  the  element.  Hances  equilibrium  cannot  be  satisfied  by  prescribing 
E  linearly  varying  axial  displacement  field  as  is  done  in  the  case  of 
the  tv’/o  noded  beam  element.  As  an  example^  a  two  noded  beam  elements 


used  to  represent  a  cantilever  loaded  at  its  free  end  with  a  concentrated 
transverse  load  Ps  yields  the  following  results  (see  Appendix  B  for 


details) : 
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These  results  are  correct  only  for  the  case  I  =1^  i.e.^ythe  reference 
axis  for  the  two  noded  beam-column  element  is  the  centroidal  axis. 

On  the  other  hand  the  same  cantilever  if  modeled  with  a  single  beam 
element  of  Fig.  l.b  yields  the  correct  results  and  ensures  a  complete 
satisfaction  of  equilibrium  in  keeping  with  the  elementary  beam  theory. 

From  this  linear  elastic  example  it  can  be  concluded  that  a  three  noded 


beam  element  is  necessary  to  analyze  the  inelastic  response  of  structures 
built  up  using  frame  elemantSs,  since  in  general  the  reference  axis  is  not 
the  neutral  axis  or  an  axis  of  constant  axial  strain  when  yielding  takes 


place  over  portions  of  the  element. 
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III.  SENSITIVITY  OF  RESPONSE  TO  APPROXIMATIONS  OF  ENERGY  ESTIMATES 


a.  Exact  Solution: 

The  exact  solution  of  a  quasi-static  inelastic  beam  problem  will 
be  developed  in  order  to  demonstrate  the  sensitivity  of  the  response  to  the 
order  of  the  integration  scheme  used  in  approximate  solutions  of  the 
problem.  Consider  the  cantilever  beam  of  rectangular  cross-section  of 
Fig.  2. a  subjected  to  a  bending  moment  at  its  tip.  The  material  of  the 
beam  is  assumed  to  be  linearly  elastic--! inearly  strain  hardening  as 
shown  in  Fig.  2.b.  For  monotonic  loading  such  a  material  is  conservative 
and  hence  the  principle  of  the  stationary  value  of  the  total  potential  energy 
can  be  used  to  obtain  the  governing  Euler-Lagrange  equations  for  the  beam. 

Euler-Bernoulli  hypotheses  imply  that 


(8) 

(9-a) 

(9-b) 

(10) 


(n-b) 


V,  is  the  volume  of  the  beam  within  which  the  beam  is  everywhere  elastic  and 

1 

is  the  volume  within  which  the  strain  at  every  point  exceeds  the  yield 
strain,  e^. 

From  purely  statical  considerations  it  is  obvious  that  the  beam  of 

2 

Fig.  2  experiences  a  constant  curvature,  — ^  and^hence,  d^  of  Fig.  3  is  a  pure 

dx"^ 

constant  and  not  a  function  of  x.  This  can  be  shown  to  be  true  by  minimizing 

the  total  potential  energy  with  respect  tosnot  only  w,but  also  d2,  assuming 

apriori  that  d2  is  a  function  of  x.  The  Euler-Lagrange  equation  resulting 

from  the  variation  of  it  with  respect  to  w  and  dp  implies  that  the  curvature, 

^-^jand  dp  are  both  constants, 
dx 

The  two  Euler-Lagrange  equations 


-d^Cj  A^dxJ-Zs  dpCj  (^)dx]-A=0 

/-I  V/  y  V 


and  the  associated  boundary  conditions 


dw(0)  _ 


tl4-a) 


(E,-E,)e  b{d?-d2)+|^[E,d^E,(d?-d^)][^^4^]=! 


2  ^2"^  3  2  2^  1  2 


(14-b) 


are  the  result  of  the  variation  of  w  with  respect  to  dp  and  w 
Equations  (13)  through  (14)  are  satisfied  by 


6 


(15) 


where 


4d,  c. 


)  Ml  -  ^)n  - 


uW  *  (fe  n  -  (^)^)] 


Use  of  Eq.  (15)  into  Eq.  (12)  yields 


ad2  =  -  / 


Substitution  of  a  from  Eq.  (16)  into  Iq.  (17)  yields  the  cubic  equation 


rp  -3(r+s)p-2(l-r)=0 


P=(^),  r=l-  ^  and  s=  -J^ - ^  (19) 

For  any  given  values  of  r  and  s,  Eq.  (18)  can  be  solved  by  trial  and  error 
to  obtain  the  corresponding  value  of  p.  The  total  strain  energy  of  deformation 
is  then  given  by 

Uj  =  [{p'^+(l-r)(l-p^)}(2ad^)‘^-3re  (l-p'^)(2ad^)-3e;r(l-p)]  (20) 


If  the  beam  is  unloaded  after  being  loaded  into  the  inelastic  range  and  if 
it  is  assumed  that  unloading  takes  place  elastically  then  it  can  be  easily 
verified  that  the  recoverable  energy  of  deformation  is  given  by 


4_3/..  s2,.2  2 


=E.,b£d^[fp"(ad^)"+r"e  "(l-p)-2ad^e  r(l-r)(l-p^)+  |(l-r)^(l-p^) (ad^)]  (21) 
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and  the  energy  dissipated  in  the  process  of  unloading  is 

W'r 

b.  Approximate  Numerical  Solutions: 

Three  numerical  solutions  to  the  problem  just  outlined  are  obtained  via 
the  minimization  of  the  total  potential  energy  using  the  three  noded  beam 
elements  of  Fig.  l.b.  and  using  numerical  integration  schemes  of  three 
different  orders  for  the  computation  of  the  strain  energy  of  deformations  [1]. 
The  first  two  schemas  use  Gaussian  quadratures  of  two  different  orders  in 
that  the  beam  is  divided  into  two  strips  in  one  case  and  into  four  strips 
in  the  other  case.  A  tvi/o  point  Gaussian  quadrature  formula  is  used  in 
each  direction  (the  length  and  the  breadth  of  the  strip).  The  third  scheme 
is  the  Newton-Cotes  scheme  which  uses  a  cubic  interpolation  with  four 
points  in  each  direction.  The  details  of  these  schemes  can  be  found  in 
any  book  on  numerical  analysis  [1].  It  must  be  emphasized  however  that 
among  integration  schemes  which  imply  that  the  integrand  can  be  approximated 
by  a  polynomial  of  some  order  the  Gaussian  quadrature  scheme  is  the  most 
efficient.  It  is  v/ell  kno'wn  for  instance  that  only  n  Gauss  points  are 
sufficient  to  integrate  a  polynomial  of  (2n-l)  degree  exactly  while  (n+1) 
points  are  necessary  to  integrate  a  polynomial  of  n-th  degree  exactly  with 
the  Mewton-Cotes  integration  schema. 

Table  1  shows  a  comparison  of  the  exact  response  with  that  predicted 
by  the  three  different  integration  schemes.  Of  main  interest  is  the  dissipative 
energy.  It  is  obvious  from  these  results  that  a  very  high  order  quadrature 
is  necessary  for  a  good  correlation  with  the  exact  solution  especially  for 
loads  where  yielding  extends  over  a  major  portion  of  the  cross-section  of 
the  beam.  The  higher  the  value  of  s  in  Table  1  the  more  extensive  the 
yielding. 
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Such  large  variations  in  the  estimates  of  deformations  and  energies  for 
a  quasi-static  case  suggest  that  a  similar  study  be  made  in  connection  with 
the  nonlinear  transient  response  of  structures.  Since  exact  solutions  to 
the  nonlinear  transient  response  of  structures  are  virtually  nonexistent  one 
has  to  rely  on  a  comparison  between  numerical  solutions  using  different 
integration  schemes. 

Figure  4  shows  the  details  of  a  solid  cross-section  beam  clamped  at  both 
ends  and  loaded  impulsively  at  the  center  by  a  triangularly  varying  impulse. 

The  material  characteristics  of  this  beam  are  assumed  to  be  as  shown  in  Fig.  2 
with  CTy  =  28.96xlO^N/m^  and  =  68. 95xlO^N/m^  and  (E2/E^ )=0.4286.  A  numerical 
analysis  of  the  transient  response  of  this  beam  is  performed  using  direct  energy 
minimization  at  each  time  step.  For  the  purposes  of  the  integration  of  strain 
energy  the  two  different  Gaussian  quadrature  schemes  shown  in  Fig.  4d  are  employed. 

Figures  5  and  6  bring  out  quite  vividly  the  differences  in  the  response 
resulting  from  the  use  of  the  two  integration  schemes.  It  must  be  remarked 
that  the  deformation- time  plots  for  the  two  schemes  are  very  nearly  identical 
and  hence  are  not  shown.  The  same  cannot  be  said  of  the  acceleration  or 
energy  plots,  however.  As  regards  the  acceleration-time  plot,  at  certain 
instances  the  accelerations  predicted  by  the  two  models  can  be  seen  to  be  off 
by  as  much  as  100%.  The  nonlinear  character  of  the  model  coupled  with  num¬ 
erical  approximations  defies  predictions  of  trends  in  response  as  evidence 
by  the  crossing  of  the  two  response  curves.  The  model  with  16  Gauss  points 
can,  in  general,  be  expected  to  provide  a  better  estimate  of  dissipation 
than  the  model  with  8  Gauss  points.  It  can  be  seen  from  Fig.  6  that  the  model 
with  16  Gauss  points  predicts  a  higher  dissipation  than  the  model  with  8 
Gauss  points  by  as  much  as  15%  over  the  interval  considered.  Again,  this 
seems  to  be  a  peculiarity  of  these  two  models  and  in  general  one  model  which 
may  at  any  given  time  predict  a  dissipation  higher  than  the  other  may  very 
likely  also  predict  a  dissipation  lower  than  the  other  at  another  time. 
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Next  the  sensitivity  of  the  response  of  thin-walled  frame  elements  is 
examined.  Such  elements  can  be  expected  to  "tone  down"  the  sensitivity  of 
response,  if  any,  to  the  order  of  the  integration  scheme.  Figure  7  shows  half 
the  finite  element  model  of  a  horizontal  box  beam  supported  on  four  uprights 
and  loaded  impulsively  at  the  center.  The  beam  and  the  supports  are  assumed  to 
be  made  from  the  same  material  the  stress-strain  curve  for  which  is  similar  to 
that  of  Fig.  2.  For  the  purposes  of  numerical  integration  using  Gaussian 
quadratures  each  wall  of  the  box  section  is  assumed  to  be  divided  into 
rectangular  strips  extending  between  the  two  nodes  of  a  frame  element  with 
six  degrees  of  freedom  at  each  node.  Table  2  lists  the  relevant  properties 
of  the  model  shown  in  Fig.  7.  Two  different  integration  schemes  are  employed  - 
one  in  which  each  wall  is  divided  into  two  strips  and  another  in  which  the 
same  is  divided  into  four  strips.  As  before  a  two  point  Gaussian  quadrature 
formula  is  used  in  each  direction.  In  the  analysis,  the  effects  of  shear 
deformations  are  ignored  in  the  interest  of  simplicity  of  the  constitutive 
model  in  the  inelastic  range.  An  impulsive  load  is  applied  at  node  3. 

Figure  8  shows  the  plots  of  acceleration  versus  time  and  Fig.  9  shows  the 
plots  of  dissipative  energy  versus  time  for  the  two  integration  schemes. 

Although  the  response  is  very  nearly  identical  in  the  initial  stages  the 
responses  for  the  two  cases  diverge  from  each  other  significantly  with  the 
passage  of  time.  As  expected  the  structure  is  much  less  sensitive  to  the 
order  of  the  integration  scheme  than  the  solid  beam  of  Fig.  4.  However,  it 
is  again  evident  that  no  trends  can  be  determined  since  both  integration  schemes 
are  only  approximate.  In  fact,  it  would  seem  that  as  a  result  of  the  approxi¬ 
mations  in  calculations  of  the  strain  energies  of  deformations  upper  bound 
solutions,  even  for  conforming  finite  element  models,  are  not  guaranteed. 

Thus  one  scheme  may  provide  a  better  answer  than  the  other  at  any  given  time 
but  the  higher  order  scheme  may  be  expected  to  give  better  results  overall. 
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IV.  CONCLUSIONS 


This  study  has  established  the  necessity  of  examining  every  finite 
element  used  for  modeling  structural  behavior  in  the  inelastic  range,  for 
consistency  of  assumptions  that  will  guarantee  satisfaction  of  equilibrium. 
Furthermore,  this  study  has  revealed  that  structures  consisting  of  frame 
members  wherein  numerical  integration  has  to  be  used  for  the  purposes  of 
evaluating  stress  resultants  or  strain  energies  will  be  sensitive  as  regards 
their  accelerations  and  energy  dissipations  to  the  order  of  the  numerical 
integration  scheme  used.  This  sensitivity  can  be  expected  to  be  more 
pronounced  for  frame  members  with  solid  cross-section  in  comoarison  with 
thin-walled  members.  This  study  demonstrates  the  need  for  higher  order 
integration  schemes  for  improved  quality  of  response  in  the  nonlinear  range. 
Such  higher  order  integration  schemes  however,  would  only  make  the  already 
expensive  numerical  analysis  of  nonlinear  response  only  more  so. 
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Table  1:  Comparison  of  Beam  Response  for  Various  Integration  Schemes 


Scheme 

Recoverable  Energy 

Dissipative  Energy 

s 

^1 

’2 

^1 

^2 

^1 

^2 

Exact 

1.122489 

1.587835 

171.19222 

269.8933 

16.0936 

105.2099 

Gauss  with 

32  stress 
ref.  pts/ 
element 

1.13313 

1.58667 

171.1086 

269.8933 

19.5555 

85.24097 

Gauss  with 

16  stress 
ref.  pts/ 
element 

1.06848 

1.63969 

170.6667 

267.6774 

0.5803 

109.20314 

Newton-cotes 
with  32 
stress  ref. 
pts. /element 

1.15733 

1.53067 

159.5500 

255.6100 

37.0460 

87.7340 

where  s  ^  =  -  0.8466,  “  1-0582 
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Table  2.  Properties  of  the  Beam  on  Elastic  Supports 


Nodes 

X 

m 

Y 

m 

Lumped  Mass  ' 
Ki log ram 

Element 

Length 

m 

^1 

cm 

^2 

cm 

t 

cm 

1 

0.8001 

0.0 

- 

1 

0.7620 

3.81 

1.463 

.0704 

2 

2.7051 

0.0 

” 

2 

0.7620 

3.81 

1.463 

.0704 

3 

0.0 

0.762 

0.07793t 

3 

0.40005 

2.9261 

7.62 

.141 

4 

0.40005 

0.762 

0.07793 

4 

0.40005 

2.9261 

7.62 

.141 

5 

1.7526 

0.762 

0.1852 

5 

0.9525 

2.9261 

7.62 

.141 

6 

0.8001 

0.762 

0.14988 

6 

0.9525 

2.9261 

7.62 

.141 

7 

'  ^ 

2.7051 

0.762 

0.11138 

tA  non-structural  mass  of  0.07793  kilogram  is  assumed  to  exist  at  node  3. 
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w(c)  =  (l-3c^+2c^)w^+£(S-2c^+C^)e^  +  (3s^-2c^)W2+J!.(-C^+C^)92 

Figure  l.a.  Two  Noded  Beam-Column  Element 


U(s)=(1-3S+25^)U^+(4c-4c^)U2+(-C+2c^)U3 
W(0  =  (l-3C^+25^)W^+£(C-2C^+S^)9^  +  (3C^-2C^)W2+)1(-C^+C^)e2 
Figure  l.b.  Three  Noded  Beam-Column  Element 
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Impulsive  load  P  (  see  Fig.  4c  belov./) 


Figure  4.  A  Solid  Cross-Section  Beam  Under  an  Impulsive  Load 
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Figure  5.  Acceleration-Time  Signature  of  Node  6  (Fig.  4a). 


Figure  6.  Dissipative  Energy  Versus  Time  (Fig.  4a) 


[o: 


Figure  7.b.  Gauss  Point  Distribution  for  a  Box  Section 


Figure  8.  Acceleration-Time  Signature  of  Node 


APPENDIX  A 

Consider  an  Eiiler-Barnoum  baam  subjected  to  enci  shears  and  inoinsnts 
For  analyzing  the  response  of  this  baarnjits  deformations  are  referenced 
with  respect  to  an  aj<is  which  does  not  coincide  with  its  centroidal  axis 
around  which  it  bands.  This  E5ns  is  designated  as  the  reference  axis. 

The  potential  energi'  expression  for  such  a  beam  is  given  by 


/  f  clf  -  z  ■^3^dAdx--Pyw(0)-P2w(A)-iv!.,||(0)-M2  ||(a)  (A-1) 


dw 


A  0 


which  Upon  performing  the  Integration  v/ith  respect  to  the  area  of  cross¬ 


section  reouces  to 


/  (-lf)^dx-EZA  /(^)(^)dx  /''(^)^dx-P,w(0)-P2w(t) 

^  0  0  dx  ^  0  dx^  ‘ 


M  dw(0) 
'  1  dx 


dw  ( £ ) 


'2  d: 


(A-2) 


where  ZA  and  I  are  respectively  the  first  and  the  second  moments  of  the 
area  about  the  reference  axis= 

Upon  requiring  that  the  variation^  otts  of  the  potential  energy  vanish 
the  necessary  Eiiler^Lagrange  equations  or  the  equations  of  equilibrium  of 
the  beam  are  obtained=  Thus^ 


t  2 

El  /  (.^::^)6(W|)dx-P^6w(0)-P^6w{£)-PU{^(0))  -  iV{47(n)=0 
0  dx^  dx^  '  2  ^  ^ 


(A-3) 
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After  integrating  by  parts  the  following  equations  and  boundary  conditions 
are  obtained. 


(A-4a) 


(A-4b) 

or 

u(0)  =  ug 
w(0)  =  wg 


^(0)-  fi* 

dx  "  ®0 

(A-5a,f) 

u(t)  =  u* 


w(£)  =  w* 


The  starred  quantities  in  the  above  equations  denote  prescribed  quantities. 
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APPENDIX  B 


a.  Cantilever  Bearn  Modeled  with  a  Tv/o-Moded  Beam  Element: 

The  two  noclecl  beam-column  element  of  Fig.  l.a.  is  used  to  model 
a  cantilever  beam  subjected  to  a  load  P  at  its  free  end.  The  potential 
energy  for  such  a  beam  Is  given  by  Eq.  (A-1)  with  P2=-P. 

It  cars  be  easily  verified  that  the  substitution  of  expressions  for  u{5) 
and  wfs)  from  Fig.  l.a.  into  this  potential  energy  expression  yields 


EAr  2  n  .2-1  EAZ  /  \/r.  A  \ 


"S'  ^^2  [1 2w^+l 2w.j0.| £-24'w.jW2+l  2w^  B2£+4e^£^- 1 2w29.j£ 


+4e.,e2£^  12w2^-12w202t+4e2£^]  +  PW2 


(B-1) 


Upon  requiring  that  the  variation  of  it  with  respect  to  the  nodal  variables 
u.|5  U2S  w.j  s  W25  and  0.js  62  be  zero  the  following  equations  are  obtained. 


Either 


or 


f^{ui-u2)-i-  ^  (e2-e^)=0 


Ul=U| 


U2=u^ 


El 


^[l  2{w.j  -W2)-!-6£(e^-e2)]=0 


w.j  =w| 
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Either 


■[12(w2-w^  )-6!!,(0^+92)]+P=O 


W2=w* 


(u2-u^)+^[6(w^-W2)+2ji(2e^+e2)]=0 


9l=e* 


(B-2a,f) 


(u^-U2)+^[6(w^-W2)+2£(9^+262)]=0 


Since,  for  the  cantilever  beam  the  nodal  displacements  u-j,  and  9-j 
are  all  prescribed  to  be  zero  the  three  equations  corresponding  to  unknown 
displacements  U2,  W2  and  02  simplify  to 


-^(l  21^2-602^^  )+P=0 


( B- 3a ,  c 


U2+^(-6w2+4e2Jl)=0 


The  above  three  equations  when  solved  simultaneously  yield  the  following 
results 


V  -Z^EI 
^  c 

W2  -  121  (3+j  ) 
c 


(B-4a,c) 


Pit 

'2  ■  2EI 
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I-AZ' 


,2,., 

Furtherrfioras  at  the  free  end  (i.e.  at  x=ft)  which  is  proportional  to  the 
dx 

bendinc  mornent  at  the  free  end  can  now  be  evaluated  to  ba 


Pi 


2EI 


(t^  -  1) 


(B-4d) 


Cl  early  5  tha  above  resul 
if  1^=1  which  is  to  say 


ts  reduce  to  the  correct  strength  of  materials  results 
the  reference  axis  coincides  with  the  centroidal  axis. 


b.  Cantilever  Beam  Modeled  with  a  Three-Moded  Beam  Element: 

The  expression  for  the  potential  energy  in  this  case  becomes 


=  lA 
21 


2,15 

1‘3 


^  Q  U 1  *  ^  U  t  U  o  “ 


2,7.2 
'2'!^ 


^  T-  "2^3 j 


^  C4LI..W, 

z 


4u2W^+3u.|e.j£~4u2e-jt+U29-]ii-4Li-,Wp+8u2W2 


-4u2Wp+u.|e2t-4u202^'^3'-i392^^ 

[  1 2w^+ 1 2w .j  e .j  £- 24w .j W2+ 1 2w -j  9 2^+  4e 1  2w2 e  1  ii+4e  .j  e 2^ ^ 

2  Z 

+Uvi'^-'l2vi^e2^  +492)i^]  +  PW2  (B-5) 

Upon  requiring  that  the  variation  of  ir  with  respect  to  the  nodal  variables 
u.j5  U2S  u^-.  W.J ,  W2  and  9.^5  02  be  zero  the  equations  corresponding  to  the 
unknown  variables  U23  u^s  1^2  ®2 
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U2-  (16u2-  IU3)  -  ■^(8w2-462JI)=0 


^3=  r  (-  I  ^3)-  ™(-4w2+3e2^)=0 


W2:  -  ™(8a2-4u3)+-^(12w2-692ii)+P=0 


2'-  -  (-4u2+3u3)+^(-6w2+4e2«.)=0 


Simultaneous  solution  of  these  equations  yields 


u.  =  -Z(§^) 


(B-6a,d) 


Z(ff-) 


(B-7a,d) 


2  ■  2EI 


These  are  identical  with  the  well  known  strength  of  materials  results 


NASA-Lansloy.  1976  CR-2732 
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